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The Cahn-Hilliard equation describes phase separation in binary liquids. Here 
we study this equation with spatially-varying sources and stirring, or advection. 
We specialize to symmetric mixtures and time-independent sources and discuss stir- 
ring strategies that homogenize the binary fluid. By measuring fluctuations of the 
composition away from its mean value, we quantify the amount of homogenization 
achievable. We find upper and lower bounds on our measure of homogenization using 
only the Cahn-Hilliard equation and the incompressibility of the advecting flow. We 
compare these theoretical bounds with numerical simulations for two model flows: 
the constant flow, and the random-phase sine flow. Using the sine flow as an exam- 
ple, we show how our bounds on composition fluctuations provide a measure of the 
effectiveness of a given stirring protocol in homogenizing a phase-separating binary 
fluid. 

PACS numbers: 47.55.-t, 64.75.+g, 47.52+j 



I. INTRODUCTION 

Phase separation and its control have received intense interest, both because of industrial 
applications [HEJE] an d the mathematics involved, in particular the Cahn-Hilliard equation. 
This equation was introduced by Calm and Hilliard [I] to describe phase separation in a 
binary alloy. Since their argument relies on the thermodynamic free energy of mixing, the 
equation is completely general, and describes any two-component system where the mixed 
state is energetically unfavorable, and where the total amount of matter is conserved. Thus 
it is used in polymer physics [2], in interfacial flows [5], and in mathematical biology [6]. In 
this paper we discuss binary liquids, and the order-parameter equation obtained by Calm 
and Hilliard describes the composition of a binary fluid. A composition c {x ) t) = indicates 
a locally well-mixed state, while a composition c{x ) t) ^ indicates a local abundance of 
one binary fluid component over another. In the context of binary fluids, it is important to 
study the influence of stirring on the two fluid components, and we therefore introduce an 
advection term into the equation. 

The time evolution of the non-advected equation has been studied extensively. In par- 
ticular, there is a proof concerning the existence and uniqueness of solutions [7j. Given an 
initial state comprising a small perturbation around the unstable, well-mixed state c = 0, 
the system forms domains of unmixed fluid that expand or coarsen in time as i 1 / 3 , the 
Lifshitz-Slyozov law [HI |9j. In many applications, the coarsening tendency of the Cahn- 
Hilliard fluid is undesirable, and this can be overcome by fluid advection, which is either 
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passive or active. In the active case p~Ql [HI [12], the composition gradients induce a back- 
reaction on the flow, while in the passive case, the composition and its gradients exert no 
effect on the flow [T3J, [HJ ESJ HE], and it is this case that we consider in the present paper. If 
the flow contains large differential shears, the coarsening is either arrested, so that domains 
grow only to a certain size [13], or can be overwhelmed entirely, so that the mixed state, 
previously unstable, is now favored [15] • This homogenization is useful in applications, for 
example in the fabrication of emulsion paints [1]. 

Stirring provides one means of controlling phase separation. Indeed, it is well known 
that a unidirectional shear flow produces banded domains, with domains aligned in the 
shear direction [HI [TTJ, [THl IE]- Other mechanisms have been employed to control the 
phase separation, including dipole interactions [20], patterned substrates [21], and surface 
forcing [22]. In [12], O Naraigh and Thiffeault study binary mixtures in thin films, and 
use the backreaction, together with surface forcing, to align the domains in a direction that 
depends on the forcing parameters. In the present paper, we focus on phase separation 
in the presence of advection and a spatially- varying persistent source. This source can 
be maintained in several ways. In [23] the spatial source is produced by thermal diffusion, 
through the Ludwig-Soret effect [24] , in which composition gradients are induced by imposed 
temperature gradients. One can also produce composition gradients simply by injecting 
matter into the system [25]. In this paper, we shall restrict our attention to the symmetric 
binary fluid, in which equal amounts of both fluid components are present. This special case 
involves a wide range of phenomena [9j EJ [19l [23], [26]. In the present work, we shall refer 
to the combination of stirring (advection) and the spatial source as 'forcing.' 

We have mentioned the often undesirable coarsening tendency of the binary fluid, and the 
efforts made to suppress it. In this paper we introduce a quantitative measure of coarsening 
suppression by studying the p th power-mean fluctuation of the composition about its average 
value. By fluctuations about the average value, we mean spatial fluctuations around the 
(constant) average spatial composition, which we then average over space and time. In 
effect, we study the time- averaged L p norm of the composition. If this quantity is small, the 
average deviation of the system about the well-mixed state c = is small, and we therefore 
use this quantity as a proxy for the level of mixedness of the fluid. An approach similar to this 
has already been taken for miscible fluids [23 [28j [29], [30], [31] . There, the equation of interest 
is the advection-diffusion equation, and fluctuations about the mean are measured by the 
variance or centred second power-mean of the fluid concentration [32], [33] . The variance is 
reduced by stirring. By specifying a source term, it is possible to state the maximum amount 
by which a given flow can reduce the variance, and hence mix the fluid. By quantifying the 
variance reduction, one can classify flows according to how effective they are at mixing. 
Just as the linearity of the advection-diffusion equation suggests the variance as a natural 
way of measuring fluctuations in the concentration, the non-linearity of the Cahn-Hilliard 
equation and its associated free energy (cubic and quartic in the composition, respectively), 
will fix our attention on the fourth power-mean of the composition fluctuations. Owing to 
Holder's inequality, a binary liquid that is well-mixed in this sense will also be well-mixed 
in the variance sense. The advantage we gain in considering the fourth power-mean is the 
derivation of explicit bounds on our measure of mixedness which have manifest flow- and 
source-dependence. 

The paper is organized as follows. In Sec. [II] we introduce the model equations and discuss 
their nondimensionalization. We introduce measures of composition fluctuations and their 
relation to fluid mixing. These measures of composition fluctuations are based on long- 
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time averages obtained from the composition of the binary fluid. Therefore, in Sec. ITTT 



we 



prove the existence of these long-time averages and find upper bounds on the measures of 
composition fluctuations. Since we are interested in minimizing composition fluctuations, in 
Sec. [IV] we obtain lower bounds on these measures. In Sec. [V] we investigate the parametric 
dependence of the upper and lower bounds for statistical, homogeneous, isotropic turbulence. 
In Sec. |VI| we compare the theoretical bounds with numerical simulations for two standard 
flows: the constant flow, and the sine flow. In the numerical simulations, we find that the 
composition fluctuations are indeed bounded by the theoretical limits we have obtained, and 
the results are dramatically different from those obtained for miscible liquids. 



II. THE MODEL EQUATIONS 

In this section, we introduce the advective Cahn-Hilliard equation with sources and 
discuss its properties. We outline the tools and notation we shall use to analyze composition 
fluctuations. For generality, the discussion takes place against the backdrop of scalar and 
vector fields in MJ 1 . 

Let c t) be the composition, that is, the scalar field c t) measures phase separation, 
with c(x,t) = indicating a locally well-mixed state, and c{x ) t) ^ indicating a local 
excess of one binary fluid component relative to the other. Let v{x ) t) be an externally 
imposed n-dimensional incompressible flow, V • v = 0, and let s (x) be a distribution of 
sources and sinks of binary fluid. The advective Cahn-Hilliard (ACH) equation describes 
the phase-separation dynamics of the scalar field c{x ) t) in the presence of flow, for prescribed 
sources and sinks, is 

dc 

— +V'\7c = DA(c 3 -c - 7Ac) +s(x). (1) 

Here D is Cahn-Hilliard diffusion coefficient and ^7 is the typical thickness of transition 
zones between phase-separated regions of the binary fluid. The finite thickness of these 
zones prevents the formation of infinite gradients in the problem. The equation is a passive 
advection equation: neither the composition nor its gradients affect the flow. 

In this paper we work with a nondimensionalization of Eq. ([T]) that leaves three control 
parameters in the problem. Therefore, we can unambiguously study limits where control 
parameters tend to zero or infinity. Let T be a timescale associated with the velocity v(x,t), 
and let Vo be the magnitude of v (x,t). Let So be the magnitude of the source variations 
and, finally, let L be a lengthscale in the problem; for example, if the problem is solved 
in a cube with periodic boundary conditions, we take the lengthscale L to be the cube 
length. It is then possible to write down Eq. ([T]) using a nondimensional time t' = t/T and 
a nondimensional spatial variable x f = x/L, 

^ + V'v • V'c = D'A' (c 3 - c - t'A'c) + S' § (a/) , 

v and s are dimensionless shape functions, Vq = TVo/L, D f = DT/L 2 , 7' = 7/ L 2 , and 
where S' Q = SqT. The quantity D' = DT/L 2 is identified with T/T D) the ratio of the 
velocity timescale to the diffusion timescale. Following standard practice, we shall now 
work with the dimensionless version of the equation, and omit the prime notation. For ease 
of notation, we shall henceforth take v to mean VqV. 

The equation (fil) has the following properties, which we shall exploit in our analysis: 
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• If the source s (x) is chosen to have spatial mean zero, then the total mass is conserved, 
d 



dt 



/ c{x,t) d n x = / .DA (c 3 - c - 7 Ac) d n x + s (x) d n x = + [boundary terms] , 



where Q is the problem domain in n dimensions and |Q| is its volume. The boundary 
terms in this equation vanish on choosing no-flux boundary conditions ft- Vc = ft- V/x = 
on dtt, or periodic boundary conditions. Here n- V is the outward normal derivative. 
In this paper we shall consider the periodic case. 



• There is a free-energy functional 

F[c]= f [t(c 2 -l) 2 + |7|V C 



8F 

d n x, /i=-— = c 3 -c - 7AC, (2) 
oc 



where \i is the chemical potential of the system [34] . For a smooth composition c (as, £), 
the free energy satisfies the evolution equation 



dF 
~dt 



D / \Vfi\ 2 d n x+ / fi(-v- Vc + s)d n x, 



and decays in time in the absence of sources and stirring. 

To study the spatial fluctuations in composition, we consider the power-means of the 
quantity c (x, t) — J Q c (x, t)d n x ) 



M p (t) = l f c ( x ^)-jk I c(x ) t)d n x 



(3) 



For a symmetric mixture in which f Q c(x, t) d n x = 0, this is simply 



M p (t) = 



\c{x ) t)\ p d n x\ = \\c 



IP? 



where we have introduced the L p norm of the composition, ||c|| p . The quantity M p is a 
measure of the magnitude of spatial fluctuations in the composition about the mean, at a 
given time. Since we are interested in the ultimate state of the system, we study the long-time 
average of composition fluctuations. We therefore focus on the power-mean fluctuations 



where (•) is the long-time average 



m p = (MV)-r.. 



(•) = lim - I (•) ds, 

t^OO t 



provided the limit exists. We shall repeatedly use the following results for the monotonicity 
of norms, 



\g(s)\ds < 



v < 11/11,, 
1 « 



t 



\g(s)\ q ds 



l<P<q, feL«(Sl), 
q>l, geL s ([0,t]) 



(4) 
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which follow from the Holder inequality. 

The Cahn-Hilliard equation and its free energy functional contain high powers of the 
composition c (c 3 and c 4 respectively), and we can therefore estimate m p for specific p- 
values. In particular, in the following sections, we shall prove the following result in n 
dimensions: 

Given a smooth solution to the ACH equation, the long-time average of the free 
energy exists, and therefore m p exists, for 1 < p < 4. 

The constraints we impose on the forcing terms are that the velocity field and its first spatial 
derivatives be bounded in the L°° norm, and that the source term be bounded in the L 2 
norm. That is, v G L°° (0, T; H l '°° (fi)) for any T G [0, oc), and s G L 2 (Q). We take our 
result one step further by explicitly evaluating upper and lower bounds for 7714, and this 
gives a way of quantifying composition fluctuations in the stirred binary fluid. 



III. EXISTENCE OF LONG-TIME AVERAGES 

In this section, we prove a result concerning the existence of the long-time average of the 
free energy, and of the power-means m p) for 1 < p < 4. 

Given the velocity field v{x,t) G L°° (0, T; H l '°° (Q)) for any T G [0,oc), the 
source s(x) G L 2 (fi), and smooth initial data for the ACH equation ([T]), the 
long-time average of the free energy exists, and thus m p exists, for 1 < p < 4. 

The proof relies on the free-energy evolution equation. Using this law, we find uniform 
bounds on the finite-time means (F) t and (Mf) t) where 

0* = 7 f(-)ds, <->=lim<.> t . 

t Jo t^oc 

Using the monotonicity of norms, the uniform boundedness of {M^) t) follows, for 1 < p < 4. 
The proof proceeds in multiple steps, which we outline below. 



Step 1: Analysis of the free-energy evolution equation 

By modifying the argument of Elliott and Zheng |[7j for the Cahn-Hilliard equation with- 
out flow and sources, it is readily seen that for smooth initial data, and for forcing terms 
with the regularity properties just mentioned, a unique smooth solution to the ACH equation 
exists, at least for finite times. Thus, we turn to the question of the long-time behaviour 
of solutions. We exploit the smoothness properties of the composition field c(x,t) and 
formulate an evolution equation for the free energy 

F H= / [\{c 2 -l) + \ 1 \Wc\ 2 ]d n x. 
Jn 

Given the smooth, finite-time solution c (#,£), we differentiate the functional F [c] with 
respect to time and obtain the relation 

dF f 90 ,n 3 

— — = / —LLCLX, U = C — C — 7AC. 
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Since c(x,t) satisfies the ACH equation ([!]), the evolution equation takes the form 

dF 



-D / |V/x|Vx+ / ii(-V'Vc + s)d n x ) 



using the no-flux or periodic boundary conditions. By averaging this equation over finite 
times, we obtain the identity 

(F) t + D( [ \Vfi\ 2 d n x\ = ( [ fisd n x\ - ( [ fiV'\/cd n x\ , (5) 
\Jn 1 1 \ Jn 1 1 \ Jn It 

We single out the quantity (F) for study. Owing to the nonnegativity of F (t), we have the 
inequality (F) > 0. Therefore, we need only consider two separate cases: (F) = 0, and 
(F) > 0. We shall show that (F) ^ is not possible, and in doing so, we shall produce a 
uniform (i-independent) upper bound on (F) t . 

Let us assume for contradiction that (F) > 0. Then, given any e in the range < e < (F), 
there is a time T E such that (F) — e < (F) t < (F) + e, for all times t > T e . Thus, for times 
t > T E) the time average (F) t is strictly positive. Henceforth, the inequality t > T £ is 
assumed. We use the condition V • v = 0, together with integration by parts, and find the 
last term in Eq. ^ becomes 

(F) t + D( [ \Vii\ 2 d n x) = ( [ fisd n x\ +( 7 / AcV'Vcd n x\ . (6) 



Now 



/ Acv-Vcd n x = - [ (d l c)[d l (v 3 d 3 c)]d n x ) 

(d l c)(d l v J )(d J c)d n x- [ (d i c)(v-V){d i c)d n x, 



'wvv w 1 1,1 



- / wWw 1 d n x, 



with 



w = Vc, Wij = \ (diVj + djVi) , 

where we have used the summation convention for repeated indices and have omitted terms 
in the integration identities that vanish as a result of our choice of boundary conditions. 
The quadratic form wWw T satisfies |tt;Vl^t(; T | < nmax^ |W#| ll^lli? which gives rise to the 
inequality 



Acv • Vcd n x 



< n ( sup | ^ | j / \Vc\ 2 d n x. 



(7) 



The matrix W is the rate-of-strain tensor. The appearance of the rate-of-strain tensor in our 
analysis shows the importance of shear and stretching in the development of the composition 
morphology. 

For each time t' G [0,t], we split the chemical potential /i into a part with mean zero, 
and a mean component: \± = Jl{t f ) + //(#,£'), where J Q ji f (as, t') d n x = 0. Then, for any 
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function (j) (sc, t) with spatial mean zero, we have the relation J Q (\>\id n x = J Q (j)n'd n x. Using 
this device, Eq. ^ becomes 



(F) t + D( / | V/xT d n x) = ( fi'sdTx ) + ( 7 / Acv • VccTa; ) . 

/ 1 \ Jn 1 1 \ Jn It 

Owing to the positivity of (F) t , we have the inequality 

D ( ( |V//| 2 d n x) <( f fjfsdTx) +(7 / Acv-VcdTx) . (8) 
\Jn 1 1 \Jn / 1 \ Jn It 

Finally, we employ the Poincare inequality for mean-zero functions on a periodic domain 

n = [o,L} 2 , 

/ T \ 2 

(9) 



Combining Eqs. 0, ([8]), and ^ gives the following inequality: 



'"\<(\M)nsh + nW t 



(10) 



where Woo = su Pt,^,i,j l^jl- There are no angle brackets around the source term because 
s (x) is independent of time. 



Step 2: Obtaining a bound on (H/i'HDt 



Using 



/ 7 I Vc| 2 d n x = [ [11' c + c 2 - c 4 ] d n x < ||// 
we obtain the inequality 

J 7 |Vc|Vx< iQ^ii/i'iyiciu + i^i^iiciii 



211CII2 + IM 



2> 



(11) 



Combining Eqs. (10) and (11) 



D 



2n 



n\2 



>*< 



/||2\ 2 

2/t 



\ s \\ 2 + n\n\* Woo<||c||^l + 2 \n\* w^iWcWi)}, 



1 

a quadratic inequality in (H/i'l^)* 2 • Hence, 



+ 2D [ 2tt ) 



/ i i \ 2 i / 2tt \ i 

[\\sh + nM-*w^(\\c\\i)?) +8D\n\-*l T ) WMi)? 



A less sharp bound is given by 



/||2\ 



^ -i^i^y ( ||s| ' 2 + n w ^ c ^y + tty w ~>^?> ^ 

which is an upper bound for (H/ZHDtj i n terms of the forcing parameters and (||c|||) t . 
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Step 3: An upper bound on 7774 



We have the free energy 



IV7 I 2 " 

Vc 


d n x= [ 




Jn 



[ [i^'OcO-Jc 4 ] d n x + \M >0. 



Hence, 



/ c 4 rf n x < 2 / c//d n x + \Q\ < 2||c|| 2 ||//|| 2 + \n\ . 
Time-averaging both sides and using the monotonicity of norms Q, we obtain the result 



|c|D t <|Q|+2|0|t(|| c ||^(||^||2)^ 



Using the bound for (H/ZIIDt m (12), this becomes 

(\\4t) t <M 



4\ 4 



\s\\ 2 + n\Q,\±W^{\\c\\\)}) +4nD\n 



1 /2tt 



Woo(\\c\\i)! 



We therefore have a ^-independent equation for the upper bound on (||c|||) t , 

{\\c\\i) t <m™(v,s,D), 
where m™ ax solves the polynomial 

(m™ ax ) 4 = I ft I 



|s|| 2 + n \fl\ 4 Wooing 



(13) 



The highest power of m™ ax on the left-hand side is (m™ ax ) 4 , while the highest power of 

3 

m max on right-hand side is (m™ ax ) 2 . Thus, this equation always has a unique positive 
solution. 

We obtain the following chain of uniform (t-independent) bounds. Each bound follows 
from the previous bounds in the chain, and the first bound follows from Eq. (13). 



• (||c|||)t is uniformly bounded, 

• (HcllDt is uniformly bounded, 

• (Ha^II!)* * s uniformly bounded, 

• (ll^lli)* i s uniformly bounded, 



(14) 
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• (F) t is uniformly bounded, 

for all t > T £ . Owing to the uniformity of these bounds, they hold in the limit t —> oc. 
The result (F) < oc implies the existence of a uniform bound for F (i), almost everywhere. 
Given the differentiability of F (i), this implies that F (t) is everywhere uniformly bounded, 
and thus, (F) = 0, which is a contradiction. Therefore, the only possibility for (F) is that it 
be zero. It is straightforward to verify that by taking (F) = 0, and making slight alterations 
in steps 1-3, the bounds in Eq. ( [Til ) still hold. 

Let us examine the significance of our result. We have shown that for sufficiently regular 
flows and source terms (specifically, v {x ) t) G L°° (0, T; if 1,00 (fi)), T G [0, oc), and s (x) G 
L 2 (^)), there is an a priori bound on the compositional free energy (F [c]). We have shown 
that the system always reaches a steady state, in the sense that (F) = 0. We have also found 
an upper bound for the measure of composition fluctuations, as the unique positive root 



of the polynomial equation Eq. (13). This bound depends only on the source amplitude, the 
diffusion constant, and the maximum rate-of-strain W^. Using the monotonicity of norms, 
this number serves also as an upper bound on m p for 1 < p < 4. Let us comment briefly on 
the volume term in the equation (m™ ax ) 4 = |f2| + .... Since this upper bound includes many 
situations, it must take into account the case where both the velocity and the source vanish. 



Then c ~ ±1 as t — ► oc, and by definition, ~ which is in agreement with Eq. (13). 



As mentioned in Sec. [1] it is desirable in many applications to suppress composition 
fluctuations, since this leads to a homogeneous mixture. In this paper, we propose advection 
as a suppression mechanism, and we would therefore like to know the maximum suppression 
achievable for a given flow. This suggests that we seek lower bounds on m p , in addition to 
the upper bounds found in this section. 



IV. LOWER BOUNDS ON THE COMPOSITION FLUCTUATIONS 

In this section we discuss the significance of the lower bound on the measure m p of 
composition fluctuations. Due to the powers of the composition that appear in the Cahn- 
Hilliard equation, it is possible to obtain an explicit lower bound for 7714, which we then use 
to discuss mechanisms to suppress composition fluctuations. After taking care of the volume 
factors, the lower bound on 7714 must be greater than or equal to the lower bound on m p) for 
1 < p < 4. Thus, a flow that suppresses composition fluctuations in the 7714 sense will also 
suppress them in in the m p sense, for 1 < p < 4. 

As discussed in Sec. [TTJ a suitable measure of composition fluctuations for a symmetric 
mixture is 

where c{x ) t) is the composition of the binary mixture and ||c|| p is its L p norm. For p — 2, 
this gives the usual variance, used in the theory of miscible fluid mixing j2?l [28j EHl I301I3T] . 
In that case, the choice p = 2 is a natural one suggested by the linearity of the advection- 
diffusion equation. In the following analysis of the ACH equation, it is possible to find an 
explicit lower bound for and we therefore use this quantity to study the suppression 
of composition fluctuations due to the imposed velocity field. Given this formula, we can 
compare the suppression achieved by a given flow with the ideal level of suppression, and 
decide on the best strategy to homogenize the binary fluid. 
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To estimate 777,4, we take the ACH equation ([!]), multiply it by an arbitrary, spatially- 
varying test function (j){x) ) and then integrate over space and time, which yields 



/ 



Dc 2 Ac() 



d n x 



sd>d n x, 



(15) 



where Q is the linear operator v • V — DA — 7DA . Using the constraint Eq. (15), the 
monotonicity of norms, and the Cauchy-Schwarz inequality, we obtain the following string 
of inequalities, 



/ 84 


cTx 


Jn 





<(||c|| 2 ||Q0 + ^A0|| 2 )<(||c||2)l 



which gives the relation 



\J a S(j)d n x\ 



We study the denominator 

-Dc 2 A<f>\\i)i < 



D(\\c 2 ml)->: 



< 



(Q<j>yd?xY +D\\A<l>U(\\c\\t)* t 



where this bound follows from the triangle and Holder inequalities. Thus we have the result 



*x\ 



or 



</ n (Q^) 2 d»x>3 +£>||A^|| 00 (||c||i> 
(Q<P) 2 d n x)" + D\\A^(\\c\\i)^ 



1 •) 

2 



> 


/ scj)d n x 




Jn 



<IMI2>* 

Using the monotonicity of norms Q , we recast this inequality as one involving only a single 
power-mean, 



(Q0) 2 cP l rr)'+£>||A0|| oo (||c|| 



> 


/ scj)d n x 




Jn 



Therefore, we have the following inequality for 777,4 - 

m 4 (q (v,D n ) + D\\A^ml)>\n\-^ [ 

Jn 



where 



q (v, D, 7) = ( / [v • V0 - DAcj) - D~fA 2 (/)] 2 d r 



"x 



X 



Thus, we obtain a lower bound for the 777,4 measure of composition fluctuations, 



m 4 > mf n , D\\ A0IU (mf n ) 3 + g («, 5, D) mf n - 



scj)d n x 



(16) 
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The cubic equation satisfied by m™ m has a unique positive root. 



To probe the asymptotic forms of (16), we rewrite the forcing terms v (as, t) and s (x) as 



an amplitude, multiplied by a dimensionless shape function. Thus, 



2\ - 



V = V V , Vb = |0| 5 (||V|| 2 

s = S s , 5 = |f2|~2 ||s|| 2j 

Then for a fixed value of So and D, and l/ > 1 (large stirring), we have q ~ 
Vb(J (£? • V0) 2 d n x)z , and the lower bound m™ m takes the form 

mmin . go IL^ n ^| y ! 

On the other hand, for fixed So and Vo, and D>1 (large diffusion), the lower bound takes 
the form 

So I ( f a §(f)d n x) I 

^ [/ n (A<£ + 7 A ! ty)Va;]*|fl|* 

It is possible to obtain similar asymptotic expressions for 777,2, by constrained minimization 
of a functional of the composition. Apart from a volume factor, the asymptotic form of 777,2 
agrees exactly with the asymptotic form of 777,4 just obtained. The functional to be minimized 
is 

$ [c] = §( / c 2 d n x^ - f (cQ(/> + Dc 3 Acf) + 50) d n x y 

where <p ( x ) is a test function. (This approach was used in [30] for the advection-diffusion 
equation.) Setting 5^ /5c = gives 



12X 2 DA6Q6 

• (17) 



6XDA4> 



Evaluation of S 2 Q/5c5c' shows that Eq. (17) produces a minimum of $[c]. Given the 
expression Q = VqV ■ V — DA — 7DA 2 , the minimum Eq. (17) is XVoV ■ V</> at large Vo- 
Substitution of this expression into the constraint (f Q cQ<p + Dc 3 A<fi + s<f> d n x) = gives 
A = - (So/V£) [(J a §(j)d n x}/(J a (v • V0) 2 d n x)] , and hence 

mmin S |(f n W| V 1S> 1 

For fixed Vo and So and large D, a similar calculation gives 

S I ( f S(f)d n x) I 



[/ n (A0 + 7 A 2 0) 2 



1 7 
2 



These expressions show that, apart from a volume factor, the lower bounds on the 777,2 and 
777,4 measures of composition fluctuations are identical in the limits of high stirring strength 
or high diffusion. In particular, the asymptotic expression 

, , 1 • So I To sSd n x\ 
m? n , Q " ' ; ' for large Vo, 
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indicates that if a flow can be found that saturates the lower bound m™, the suppression 
of composition fluctuations can be enhanced by a factor of V^ 1 at large stirring amplitudes. 
Such a flow would then be an efficient way of mixing the binary fluid. 



V. SCALING LAWS FOR m 4 



In this section we investigate the dependence of the measure of composition fluctu- 
ations on the parameters of the problem, namely the stirring velocity v ) the source s, and 
the diffusion constant D. For simplicity, we shall restrict our interest to a certain class of 
flows, which enables us to compute long-time averages explicitly. 

The lower bound for the measure of composition fluctuations is the unique positive 
root of the polynomial 



s6d n x 



= 0, 



(18) 



where (j) (x) is a test function and 

q 0, D,-y) = (J [v • V0 - DA0 - D 7 A 2 0] 2 d n x 

Following [271 EE], we specialize to velocity fields whose time average has the following 
properties, 

K 2 



(vi (x, •)) = 0, (vi (x, •) vj (x, •)) = —Sij. 



n 



(19) 



The flow v (x,t) is defined on the n-torus [0, L] n . A statistically homogeneous and isotropic 
turbulent velocity field automatically satisfies the relations (19), although it is not necessary 



for v to be of this type. The source we consider is monochromatic (that is, it contains 
contains a single spatial scale) and varies in a single direction, 

s = V2S sin (k s x) . (20) 
Our choice of velocity field makes the evaluation of q (v, D,^) particularly easy: 



q (v,D,7) 



V 2 
n 



\V<P\ 2 d n x + D 2 / (A0 + 7 A 2 ^ 
[o,L] n AoM n 



d n x 



In studies of the advection-diffusion equation [30], it is possible to find an explicit test 
function (f) that sharpens the lower bound on m 2 . The procedure for doing this depends on 
the linearity of the equation. Here, this is not possible, and for simplicity we set (j) = s. This 
choice of cp certainly gives a lower bound for the 1714 measure of composition fluctuations, 
with the added advantage of enabling explicit computations. Having specified the coefficients 



of the polynomial in Eq. (18) completely, we extract the positive root of this equation, and 
find the lower bound m™ ax , as a function of Vq. The results of this procedure are shown in 
Fig.[T} 
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FIG. 1: (a) The lower bound for mf m as a function of Vo for a monochromatic source; (b) the 
dependence of mj 1111 on the velocity amplitude Vo. In (a), the scale of the source variation decreases 
in integer multiples from A: s = 2tv/L in the uppermost curve, to k s = 8tv/L in the lowermost curve, 
while in (b) the source scale is set to 2tv/L. In both figures, we have set D = Sq = 1. 



Let us examine briefly the scaling of the upper bound m 
The upper bound satisfies the polynomial equation 



max 
4 



with the problem parameters. 



(m 4 max ) 4 = |fi| 



2\Q\ 



D 



2tt J 



So\Q\z +n\Q 



\ 2 , ,i /27T 

W 00 m 4 max j +4nD\Q\* I — 



Woo (m 4 max ) 2 



(21) 



which depends only on the diffusion D, the source amplitude So, and the maximum rate- 

i 

of-strain W^. For large, the flow-dependence of the upper bound is m™ ax ~ W£>. This 
dependence is verified by obtaining the positive root of Eq. (21), which is a function of W^. 
The results are shown in Fig. [2| 

In this section we have investigated the parametric dependence of the theoretical bounds 



max, mm 



for flows with the properties in Eq. (19). We note that for any nonzero source, 



the lower bound mf m is nonzero, meaning that no matter how hard one stirs, there will 
always be some inhomogeneity in the fluid, and this is in fact true for any flow. However, 
the number m™ m tells us how much homogeneity we can achieve and is therefore a yardstick 
for stirring protocols. We use this yardstick to test model flows in the next section. 



VI. NUMERICAL SIMULATIONS 



In this section we solve Eq. ([!]) numerically for two flows, and verify the bounds obtained 
in Sees. III-V We use the sinusoidal source term in Eq. (20) with periodic boundary condi- 
tions, and the source scale k s therefore takes the form (2ir/L) j, where L is the box size and 
j is an integer. We specialize to two dimensions and study two standard flows that are used 
in the analysis of mixing: the random-phase sine flow [25j [29j |35j |36l |37j, and the constant 
flow [31]. 
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FIG. 2: The dependence of the upper bound mj 1 ^ on the maximum rate-of-strain Woo- The source 



affects the upper bound only through its square mean So = \Q\ 



s 2- 



Random-phase sine flow 



The random-phase sine flow is the time-dependent two-dimensional flow 



(x, y,t) = y/2 V sin (k v y + fa) 



jr<t< (j + \) t. 



(x, y,t) = V2 V sin (fc v x + ^) , v x = 0, ( j + |) r < t < (j + 1) r, 



(22) 



where 0j and ^ are phases that are randomized once during each flow period t, and where 
the integer j labels the period. The flow is defined on the two-dimensional torus [0, L] . 
The time average of this velocity field has the properties listed in Eq. (19). Because of its 



simplicity, the sine flow is a popular testbed for studying chaotic mixing [25, 29, 1351 1361137] . 

We solve Eq. ([j]) with the flow in Eq. (22) using an operator splitting: the advection 
step is carried out using the lattice method of Pierrehumbert j25j [37], and the subsequent 
Cahn-Hilliard and source steps are implemented simultaneously using a spectral method [9]. 
The nondimensionalization outlined in Sec. [TT] is appropriate here: the unit of time T is 
identified with the is the flow period r, and the unit of length is the box size L. The control 
parameters in the problem are the dimensionless velocity Vo, the dimensionless diffusion D, 
and the dimensionless source amplitude So = |fi| _ 2||s|| 2 . We use Vo as a measure of stirring 
intensity and fix the other parameters in what follows. The flow we choose is chaotic at all 
stirring amplitudes and given our choice of scaling, has Lyapunov exponent 



A- 0.2361/o 2 , Vo < 1; A - log 



V 2 



The lattice method with its splitting of the advection and diffusion steps, is effective only 
when diffusion is slow compared to advection, that is, T/T D <C 1. We therefore set D = 10 -5 , 
with r — L = 1. A numerical experiment with Vo = shows that So = 5 x 10 -4 gives rise 
to a morphology that is qualitatively different from the sourceless case, and we therefore 
work with this source amplitude. Finally, following standard practice [HI H3] . we choose 
7 rsj Ax 2 , the gridsize. 

Using these new scaling rules, and the identity = y2Vok Y for the sine flow, we recall 



15 



the large-stirring forms of m™ ax ' mm , For Vq ^> D, the lower bound has the form 

So J a sin (k s x) <j)d 2 x 



1 •) 

2 



(23) 



with the power-law relationship m™ m ~ V 1 , while for small Vo 3> 1 the upper bound has 
the form 

i 

V^. (24) 



TT 2 D 

These scaling results are identical to those for the advection-diffusion problem |29|. 

Before studying the case with flow, we integrate Eq. ([I]) without flow, to verify the effect 
of the source. For a sufficiently large source amplitude, the composition phase- separates 




(a) 




(b) 



(c) 



(d) 



1.5r 



0.5 



1 2 3 4 

10 10 10 10 

t 

(e) 

FIG. 3: The composition of the binary fluid for S = 5 x 1CT 4 and (a) t = 100; (b) t = 1000; (c) 
t = 2000; (d) t = 8000. A steady state is reached in (d), evidenced by the time dependence of 7774 
in (e), where 7714 is constant for t ^> 1. 

and forms domains rich in either binary fluid component. These domains are aligned with 
variations in the source. A steady-state is reached and 7714 attains a constant value, as 
seen in Fig. |3| On the other hand, for small source amplitudes, we have verified that that 
the domains do not align with the source, and their growth does not saturate. We do not 
consider this case here, since we are interested in sources that qualitatively alter the phase 
separation. These different regimes are discussed in [23] . 

We consider the case with flow by varying Vo, and find results that are similar to those 
found in [15], for the same stirring mechanism without sources. For all values of Vq, the 
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composition reaches a steady state, in which ||c||4, the pre-averaged form of 7714, fluctuates 
around a constant value. For small values of Vb, the domain growth is arrested due to 
a balance between the advection and phase-separation terms in the equation, while for 
moderate values of Vb, the domains are broken up and a mixed state is obtained. At large 
values of Vb, the 777,4 measure of composition fluctuations saturates: further increases in Vb 
do not produce further decreases in 777,4. At these large values of Vb, the source structure is 
visible in snapshots of the composition, as evidenced by Fig. [4} 

We investigate the dependence of composition fluctuations on the stirring strength Vb, 
and show the results in Fig. [5| The theoretical upper and lower bounds on 777,4 depend on 
Vb and are obtained as as roots of Eqs. (18) and (21). In the limit of large Vb, these bounds 

~L 

have power-law behaviour, with m^ ax ~ V 2 and m™ m ~ Vb" 1 , as demonstrated by Eqs. (23) 



and (24). The numerical values of 777,4 are indeed bounded by these limiting values, although 
the Vb-dependence is not a power law. Instead, the function 777,4 (Vb) is a non-increasing 




J 



(a) 



(b) 



(c) 



FIG. 4: A snapshot of the steady-state composition for (a) Vb = 0.001; (b) Vb = 0.1; (c) Vb = 10. 
In (a) domain growth is arrested, in (b) the domains are destroyed and the binary fluid mixes, while 
in (c) 7774 measure of composition fluctuations is minimized, and the source structure is visible. 



function, with a sharp drop occurring in a small range of Vb-values. Thus, the fluid becomes 
more homogeneous with increasing Vq. We discuss the effect of stirring on the inhomogeneity 
of the fluid by introducing the notion of mixing enhancement. 

We measure the ability of a given stirring protocol to suppress composition fluctuations 
by the mixing enhancement. Similar ideas are often applied to the advection-diffusion equa- 
tion [271 EHl EHJ [30j EI] . We define the dimensionless mixing enhancement 



mf n (Vb = 0) 



_ p 



rn, 



(Vo) 



For a given flow, the number r\ v quantifies the flow's ability to suppress composition fluctua- 
tions. In a well-mixed flow, the local deviation of c (x : t) away from the mean will be small; 
a small m^-value is a signature of a well-mixed flow. We are therefore justified in calling 
r] p the mixing enhancement. We obtain some control over the mixing enhancement 774 from 
the inequalities of Sees. [TTT]and[lVj Based on these inequalities, the mixing enhancement is 
bounded above and below, 



111111 

74 — 



mf n (Vo = 0) 
m^ ax (Vb) 



<V4< 



mf n (V = 0) 
mf n (Vo) 
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FIG. 5: (a) The 7714 measure of composition fluctuations or mixing for the sine flow, as a function 
of the stirring parameter Vq. The values of D and Sq are given in the text. The upper and lower 
bounds are shown for comparison; (b) The mixing enhancement 774 for the sine flow, with the upper 
and lower bounds shown for comparison. 



We have plotted the upper and lower bounds on the mixing enhancement for the case of 
monochromatic sources in Fig. [5](b). The maximum enhancement always exceeds unity in 
this case, which implies the possibility of finding stirring protocols that homogenize the 
fluid. On the other hand, the minimum enhancement is less than unity, which indicates the 
possibility of finding stirring protocols that actually amplify composition fluctuations, and 
this amplification depends weakly on the maximum rate-of-strain W^. This latter case is 
not surprising, given that a uniform shear flow causes the domains of the Cahn-Hilliard fluid 
to align, rather than to break up. The sine-flow enhancement is a non-decreasing function 
of the stirring parameter Vq. At small values of Vb, small increases in the vigor of stirring 
lead to small small increases in the mixing enhancement. There is a window of intermediate 
Vo-values for which the mixing enhancement increases sharply with increasing Vo- At higher 
values of Vb, the efficiency saturates, so that further increases in the vigor of stirring have 
no effect on composition fluctuations. The saturation is due to finite-size effects: the sine 
flow wraps filaments of fluid around the torus as in Fig. He). 



Constant flow 



We study the flow (v x , v y ) = (0, Vb), where Vb is a constant. We choose a nondimension- 
alization that is set by the diffusion time = L 2 / and obtain the following parametric 
version of Eq. Q, 

% + = A' (c 3 - c - YA'c) + V2S> sin (k'x') , (25) 

where Vq = LVq/D, 7' = 7/L 2 , and So = SqL 2 /D. We immediately drop the primes from 



Eq. (25). We fix 7 and So and vary the flow strength Vb- This flow does not satisfy the 
time-correlation relations ( fl9| ), although the maximum rate-of-strain has the simple form 
Woo = 0. Using this information, and the test function (j) = s(x), the upper and lower 
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(a) (b) (c) 

FIG. 6: A snapshot of the steady-state composition for (a) Vq = 10; (b) Vq = 100; (c) Vq = 1000. 



bounds obtained in Eqs. (13) and (16) are independent of the flow strength. 



We solve Eq. (25) numerically for various values of Vq and present the results in Fig. \6\ 



For small Vb, the composition morphology is similar to the flowless case seen in Fig. [3j( 
except now the domains are uniformly advected in a direction perpendicular to the source 
variation. The small- Vb case is shown in Fig. |6](a). As Vb increases, the domain boundaries 
are distorted due to the advection, while for large Vb, the advection is sufficiently strong to 
break up the laminar domains. The domain-like structure persists at late times however, 
and narrower laminar domains form. 

The m4 measure of mixedness is almost constant across the range of stirring parameters 
< Vb < 1000. For Vb = 1000, is slightly smaller than its value at Vb = 0, due to 
the presence of more interfaces. This difference is small however, and increasing Vb does 
little to mix the fluid. This is not surprising, given that local shears are necessary to break 
up domain structures [llj, and that such shears are absent from constant flows. What this 
example shows however, is the difference between a miscible mixture, and a phase-separating 
mixture. For a diffusive mixture with the sinusoidal source we have studied, the constant 
flow discussed here is optimal for mixing [3]]; for a phase-separating mixture, the constant 
flow badly fails to homogenize the mixture. 



VII. CONCLUSIONS 

We have introduced the advective Cahn-Hilliard equation with a mean-zero driving term 
as a way of describing a stirred, phase-separating fluid, in the presence of sources and sinks. 
By specializing to symmetric mixtures, we have studied a more tractable problem, although 
one with many applications. 

Our goal was to investigate stirring protocols numerically and analytically, and to deter- 
mine the best way to break up the domains in the Cahn-Hilliard fluid and achieve homog- 
enization. To this end, we introduced the m p measure of composition fluctuations. Since 
in a well-mixed fluid, the composition exhibits spatial fluctuations about the mean, with 
better mixing leading to smaller fluctuations, we used m p as a measure of mixedness or 
homogeneity. We proved the existence of m p for long times, for 1 < p < 4, and obtained a 
priori upper and lower bounds on 7714, as an explicit function of the imposed flow v{x ) t) ) 
and the source s (x). 

We compared the level homogeneity achieved by the random-phase sine flow and the 
constant flow with the lower bound, and found that the sine flow is effective at homogenizing 
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the binary fluid, while the constant flow fails in this task. This is not surprising, since it is 
known that differential shears are needed to break up binary fluid domains, although it is 
radically different from the advection-diffusion case, where the constant flow was the optimal 
mixer. The question of whether or not a flow is a good mixer in this context was discussed 
using the mixing enhancement, defined in Sec. [VTJ Given such a definition, it is possible 
to compare stirring protocols and find the optimal protocol for mixing a binary fluid. Our 
upper bound on the enhancement provides a meaningful notion of this optimality. This 
result may be useful in applications where the homogenization of a binary fluid is desirable, 
since we have set a lower limit on precisely how much homogeneity can be achieved. 
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